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Abstract 

Let c m , with m = 1,..., N (with N an arbitrary positive integer, N > 2) be 
the N zeros (arbitrarily ordered) of the Hermite polynomial iJjv (c), of order N 
and argument c: Hjy (c m ) = 0. Let the monic polynomial pjy (z) of degree N 
in the variable z be defined as follows: 

N N 

PN (z) =/ + ^ (c m Z N ~ m ) = JJ (z - Z n ) ■ 

m =1 n= 1 

The first equality identifies the N coefficients c m of this polynomial p^ ( z ) as 
the N zeros c m of the Hermite polynomial of order N; note that there are 
N\ such different polynomials pn ( 2 ), depending on the ordering assignment 
of the N zeros c m of the Hermite polynomial of order N. The second equality 
identifies (uniquely up to permutations) the N zeros z n of the polynomial pjv (z). 
In this paper we define in terms of these N zeros z n two N x N matrices 
M* 1 ) and Ml 2 ) with integer respectively square-integer eigenvalues Ah'* = Tn 
respectively = m 2 , m = 1, The technique whereby these findings are 

demonstrated can be extended to other named polynomials. 

Keywords: zeros of polynomials; Hermite polynomials; special functions, 
Diophantine matrices. 

MSC: 11C08, 70F10, 70K42, 11D41, 33E99. 

1 Introduction 

Notation 1.1. Unless otherwise indicated, hereafter N is an arbitrary positive 
integer , N > 2, indices such as n, to, £, j, k,... run over the integers from 1 to N, 
boldface lower-case letters indicate IV-vectors (for instance the vector z has the 
N components z n ) and boldface upper-case letters indicate N x N matrices 
(for instance the matrix M has the N 2 components M nm ). The imaginary 
unit is hereafter denoted as i (i 2 = — 1; of course i is not a N- vector!). For 
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quantities depending on the independent variable t a superimposed dot indicates 
differentiation with respect to t. The Kronecker symbol S nm has the usual 
meaning: S nm = 1 if n = to, S nm = 0 if n ^ to. And we adopt throughout the 
usual convention according to which a void sum vanishes and a void product 
equals unity: J2f=j fj = 0’ TIy=j fj = 1 if J > AT. Finally we introduce the 
following two convenient notations: 

M z ) = E z Sl z S2 ---z Sj , (la) 

l<S 1 <S2<...<S j <N 

°rn,j{ Z) = 51 (lb) 

where of course the symbol 


E 

l<Sl<S2<...<Sj<N 

denotes the sum from 1 to AT over the j integer indices Si, S 2 , ■ ■ ■, Sj with the 
restriction that si < S 2 < ... < Sj, while the symbol 

E 

l<s 1 <s 2 <...<s j - 1 <N ; Sfc^ra, k=l,...,j—l 

denotes the sum from 1 to AT over the (j — 1) indices si, S 2 , ■ ■ ■, Sj -1 with the 
restriction that si < S 2 < ... < Sj -1 and moreover the requirement that all 
these indices he different from m. We note that u m| i(z) = 0 according to the 
convention (see above) that a sum over an empty set of indices equals zero. ■ 
A class of properties satisfied by the N zeros z n of several named polyno¬ 
mials of order N has been obtained via the identification of N x N matrices, 
constructed with these N zeros z n , which feature eigenvalues displaying re¬ 
markable Diophantine properties: see for instance nm. a general technique 
to arrive at such results goes through the following two steps. 

First step. A dynamical system is manufactured, characterized, say, by N 
Newtonian equations of motion (written as follows in AFvector notation), 

C = i f (C) i C = CW, (2a) 

featuring the following two peculiar properties, (i) This dynamical system is 
isochronous with period 2i r: all its solutions are completely periodic with period 
2 tt, 

c (t + 2 tt) = C (t) . (2b) 

(ii) This dynamical system features the equilibrium configuration 

C(t)= z , (3) 

where the N components z n of the A"-vector z are just the N zeros of the named 
polynomial under consideration. 
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Second step. The behavior of dynamical system (I2al) in the infinitesimal 
neighborhood of its equilibrium (0) is investigated in the standard manner, by 
setting 

C(f) = z + ev(t) (4) 

with e infinitesimal, and by thereby obtaining from Cal) the linearized system 


v(t) = i M(z) v(f) , 


(5a) 


where of course the (f-independent) TV x TV matrix M (z) is defined componen¬ 
twise as follows: 

(5b) 
C=z 


M nm (z) = 


d fn (C) 


9 Cv 


Hence the solution of this linearized system © is a linear combination, with 
constant coefficients, of exponential functions exp (iA m f), where A m are the TV 
eigenvalues of the TV x TV matrix M (z). But if the f-evolution of the TV solutions 
(t) of the dynamical system Cal) are characterized by the periodicity property 
(l2bl) . the TV quantities v n (t) —see © must also possess this same property. 
Hence the TV eigenvalues A m must feature the Diophantine property to be inte¬ 
gers ; and these integers can be actually identified by comparing the behavior of 
the solvable system with that of its linearized version in the infinitesimal and 
immediate vicinity of the relevant equilibria. 

Remark 1.1. An analogous process can be applied to a second-order dy¬ 
namical system C, = f(C) satisfying conditions (i) and (ii), instead of the first- 
order system (I2al) . In this case, the linearization obtained in the second step 
leads to a second-order linear system v(f) = —M(z)v(f), where the matrix M(z) 


C=z 


is defined componentwise by M nm (z) = — d J^' > 

In this paper we again follow this approach, but with a new twist—based on 
a new development allowing to manufacture larger classes of solvable dynamical 
systems 013 — yielding results which seem new and indeed somewhat surprising. 
Indeed, the TV x TV matrix (z), which we identify as having TV integer 
eigenvalues, is defined in terms of the TV zeros z n of (monic) polynomials of 
degree TV in z the TV coefficients of which are the TV zeros of Hermite polynomials 
of degree TV. Moreover, by starting from the system of Newtonian (second- 
order) ODEs the solvability of which has been demonstrated in [5] rather than 
from a first-order system of ODEs such as Cal) , see Remark 1.1, we provide 
an analogous derivation of another TV x TV matrix M^ 2 ) (z)—again constructed 
with the same zeros z n —featuring TV squared-integers as eigenvalues. These 
findings are reported in the following Section 2 and proven in Section 3. Section 
4 (“Outlook”) outlines tersely some possible further developments. 


2 Results 

In this section we report our findings; their proof is provided in the following 
Section 3. 


3 






Proposition 2.1. Let c m with m = 1 be the N zeros (arbitrarily 

ordered) of the Hermite polynomial Hn (c) of order N and argument c (see for 
instance 0): 


H n (c) = N\ 


[[AT/2]] 

E 


k=0 


(-l) k ( 2 c) N ~ 2 
k\ (N — 2 k)\ 


k 


(6a) 


H n (c m ) = 0 , m=l,...,N , (6b) 

where of course the notation [[AT/2]] denotes the integer part of N/ 2, i. e. N/2 
if N is even , {N — 1) /2 if TV is odd. Let the rnonic polynomial p^ (z) of degree 
N in the variable z be defined as follows: 

N N 

PN (z)= z N +j 2 ( c - zN ~ m ) = n (* - • ( ? ) 

m=1 n =1 


The first equality identifies the N coefficients c m of this polynomial p n (c) as 
the N zeros c m of the Hermite polynomial of order N ; note that there are N\ 
such different polynomials pn (z) = p(z), pi = 1 depending on the 

ordering assignment of the N zeros c m of the Hermite polynomial of order N. 
The second equality identifies (uniquely up to permutations) the N zeros z^ of 
the polynomial p^ ( z ). 

Let the N x N matrices = MA 1 ) (z^) be defined componentwise 

as follows in terms of the N zeros (z[^\ ..., z^ ) = z^ of the polynomial 


P { N 0) : 


( z w) 


fr _I_ 

II r>) _ >), 

. Y~l J . 


f r 

N-j 

[ 4 »>] 



+ 

3,m ^ 


N 

E 

s=l,s^j 


„(m) 


» 


, O) 

Wj m - W s,m 

EFp 


(8a) 


where 


and 



(8b) 

(8c) 


(see Notation 1.1). 

Then the N eigenvalues A^ of the N x N matrices M^^are given by the 
following neat ( Diophantine) formula: 


A^ )= m, m= 1,2,..., TV . ■ 


(9) 


Remark 2.1. Note that in Proposition 2.1 all the different N\ matrices 
MTXaO = M (1 ) (z^)), where p = 1,2,..., TV!, feature the same set of N eigen¬ 
values A^. Also note that, although definition (15al) of the matrix de¬ 

pends on the ordering assignment of the N zeros z^ of the polynomial p^\z), 
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such different assignments produce the same matrix up to a reshuffling 

of its lines and columns and have therefore no relevance for the eigenvalues of 
M< : indeed, a switch of the two zeros and z^ of the polynomial p^ ( z ) 
yields a new matrix which corresponds to the original matrix 

up to the exchange of the n-th and the m-th rows and columns of (see 

also below Remark 3.2). ■ 

Proposition 2.2. Assume the notation of Proposition 2.1. If the N x N 

matrices M/ 2 )^) = M/ 2 )(z m ) are defined by 


= MW (zW) 


N 

•E 

3 =1 


An) 


lN-j 


» 


N 

TT 

1 1 

11 

[Ar 1 - J 

N 

E 

E 

_i 

V 

UJjm 'Ws,m 


1 An) (a0u 

s=l,sj£j \~j L ' s J . 

:s are given by 

2,..., 

N. 


( 10 ) 


( 11 ) 


Remark 2.2. Note that in Proposition 2.2 all the different N\ matrices 
= M^ 2 ) (z^)), where p = 1, 2,.. ., TV!, feature the same set of N eigen¬ 
values A^). Also note that, although definition (flT)l) of the matrix M^ 2 )^) de¬ 
pends on the ordering assignment of the N zeros zif' 1 of the polynomial p^\z), 
such different assignments produce the same matrix up to a reshuffling 

of its lines and columns and have therefore no relevance for the eigenvalues of 
Ml 2 )!^) (cf. Remarks 2.1 and 3.2). ■ 

Examples 2.1 and 2.2 illustrate Proposition 2.1 for the cases where 
N = 2 and TV = 3. 

Example 2.1. Let us construct the 2! = 2 matrices (see Propo¬ 

sition 2.1) for the case N = 2, where the index p £ {1,2} identifies the 2 
ordering assignments of the 2 zeros ^~^ j of the Her- 

mite polynomial of degree 2, ^(c) = 4c 2 — 2. In this special case the 2 numbers 
c^\ enter in definition (18al) of the matrix only via the expression 

[c^ — 4 m) ]~ 2 = 1/2, which does not depend on the index p distinguishing the 
two different orderings of the two zeros cf 1 ' 1 of the Hermite polynomial 

H 2 (c). Therefore, both choices of the pair yield the same formula 

for the matrix MWW in terms of the zeros of the polynomial 

z 2 + z + : 




1 - 2 . 


(e) 


2[zf' > —z. 

i-i 4 m) ] 2 / 


(m) 

2_ 

(p)i 


2 


W ] 2 
(m)i 




(12a) 


5 















( 12 b) 


Here and 


are of course the 2 zeros of the polynomial 


pM(z) = z 2 + 


(-1 r (i—~) 

V2 


p = 1 , 2 ! = 1,2 


hence 


z ± 


(-If ± yj 1- (-If 4 V2 


2f2 


M = 1,2 


( 12 c) 


(note that for p = 1 these two zeros are real , for /r = 2 they are complex 
conjugate). 

It can be easily verified that each of the matrices Mf'W, ^ = 1 , 2, given by 
formula (11 2 all has the two eigenvalues 1 and 2 consistently with Proposition 
2 .1. And it is easy to see from formula ( 112 al) that the matrix Mf) (z^, z[^) can 
be obtained from the matrix Mf f zf"\ z^) by switching its rows and columns, 
consistently with Remark 2.1. ■ 

Example 2.2. For N = 3 the 3 zeros of the Hermite polynomial H 3 (c) = 
8 c 3 — 12c = 4c (2c 2 — 3) are 0, ± y/3/2. So we make the following 6 = 3! different 
assignments: 


,(1) _ 
1 ■ 

= 0 , <4^ 

= \/372 , 4 1) = 

-y/ 3 / 2 , 

(13a) 

,( 2 ) _ 
1 “ 

= 0 , c? 

= -Vw, 4 2) 

= x/3/2, 

(13b) 

,(3) _ 
1 “ 

= VW, 

4 3) = 0 , 4 3) = 


(13c) 

,(4) _ 
1 “ 


4 4) = -V 372 , 

4 4) = 0 , 

(13d) 

,(5) _ 
1 " 

= -y/3/2, 

4 5) = vWz, 

0 

II 

iSco 

0 

(13e) 

.(6) _ 
1 " 

= -x/372 , 

4 6) = 0 , 4 6) 

= 

(13f) 


and correspondingly we define the following 6 polynomials p 3 ^ (2), of third 
degree in z, and their sets of 3 zeros z[^\ z^\ z 3 ^: 


P 3 


(m) 


(f = z 3 + [ ( 


.00 ^3-r 


n [-- - 


jt = 1, ...,6 . (14) 


The 3 zeros z[^\ z% , z^ can be easily obtained from this formula (if need 
be, via Mathematica), but rather than writing here their exact values (involv¬ 
ing square and cubic roots) we simply provide their (of course approximate) 
numerical values in decimal form: 

z[ 1] = 0.7090, z^ = -0.3545-1.2656 i , z^ = -0.3545 + 1.2656 i , (15a) 

z[ 2) = 0.7202 - 0.5758 i , = 0.7202 + 0.5758 i , zf 1 = -1.4405 , (15b) 

z{ 3) = -1.0031 + 0.7492 i , ^ 3) = —4.0031 — 0.7492 i , z^ = °- 7814 , ( 15c ) 
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zf ] = 0 , 4 4) = —1-8772 , 4 4) = 0.6524 , (15d) 

z[ 5) = 0 , 4 5) = -1-8772 , 4 5) = 0.6524 , (15e) 

z[ 6) = -0.7814 , 4 6) = 1.0031 - 0.7492 i , z^ = 1.0031 + 0.7492 i . (15f) 

Of course these numbers are defined up to permutations, but, as explained 
below, this has no relevance for the eigenvalues of the matrices and 
see Remarks 2.1 and 3.2. Note moreover that, as it happens, the zeros of the 
polynomials p^\z) and p^\z) are the same. 

From formulas (I 8 cl) for N = 3 we moreover obtain that 


„(/d 


An) 


An) 


„(m) 


_ _1 .An) _ An) i Am) .An) _ _ in) O) _ i o o j 

“'l ,m ~ B “ l 2,m “ "m+l tz m+2i w 3,m — z m+l z m+2 > _ -L z , ° muu 1°) ■ 

( 16 ) 

By inserting these values of the quantities into formula (15al) for iV = 3 
we obtain the components of the corresponding 6 matrices 

=- 1 

- 1 V - K - r) 


rj/d » 


An) An) 

Zn - Z n+1 


Z n z n+2 


( 


2 

1 _ 

An) 

^ n 

+ z (n) 

' n 

l 




(fO 


An) 


ra +1 ‘ X/ m +2 


An) An) 

c m -\-1 ^m+2 




- 1 



( 

An) 

~n 

i _i_ An) i An) 

1 ' /C m+ 1 ' ^m+ 2 


\ 


LM m] 
C 1 c 2 

2 


+■ 


_i -l An) 2 W 

1 i ^ra +2 


»(#*) 


+■ 


Am) 

■'m +1 


Am) _ (m) 
6-1 


_(m) , Am) _(m) 

4 m+2 ‘ ^ra+l ^m+2 


c (m) _ c (m) 



m = l,2,3 mod(3) . (17) 


It can be checked by direct computation that the 6 matrices given by m 
all feature the 3 eigenvalues 1, 2 , 3. And it is again plain that a permutation 
of the zeros of the polynomial p^\z) results in a corresponding permutation of 
the rows and the columns of the matrix M^Am), see Remarks 2.1 and 3.2. ■ 
Let us end this section by displaying the matrices M^ 2 ) (z), defined in Propo¬ 
sition 2 . 2 , for N = 2 and N = 3. Note that these matrices M^ 2 ) (z) do not 
coincide with the matrices [MA) (z)] (see those reported in Examples 2.1 
and 2.2), although they of course feature the same squared-integer eigenvalues. 

Example 2.3. Let us construct the 2! = 2 matrices M( 2 )(m) f 01 - the case 
N = 2, where, as in Example 2.1, the index p £ {1, 2} identifies the 2 ordering 

assignments of the 2 zeros (c^, ''j of the Hermite poly¬ 

nomial of degree 2, # 2 ( 0 ) = 4c 2 — 2. Similarly to Example 2.1, in this special 
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case the 2 numbers c^\ enter in definition (flUl) of the matrix M/ 2 )^) only 
via the expression (c£^ — C 2 ) -4 = 1/4, which does not depend on the index 
fi distinguishing the two different orderings of the two zeros c^\ of the 
Hermite polynomial /^(c). Therefore both choices of the pair (c [ ll> , yield 

the same formula for the matrix in terms the zeros = {z^, of 

the polynomial z 2 + z + : 


M (2 ^) = 


3 1-2 


£2 s 1 -^) 2 \ 


i-(2^y 


5 , 3 i~g£2 4*° 

2 ' 2 -,(/*)_-(/*) / 
Z 1 z 2 / 


(18) 


Here and z^ are of course the 2 zeros of the polynomial (I12bl) given by 
(112cl) (note, again, that for /r = 1 these two zeros are real , for /r = 2 they are 
complex conjugate). 

It can be easily verified that each of the matrices M/ 2 ^), fi = 1, 2, given by 
formula (fTSl) has the two eigenvalues 1 and 4, consistently with Proposition 
2 . 2 . And it is easy to see from formula (1T5T) that the matrix M^ 2 )(^ ^U^can 
be obtained from the matrix M/ 2 )(z^, z^) by switching its rows and columns, 
consistently with Remark 2.2. ■ 

Example 2.4. For N = 3 the 3 zeros of the Hermite polynomial (c) = 
8 c 3 — 12c = 4c (2c 2 — 3) are 0, ± y/3/2. As in Example 2.2, we make 3! = 6 
different assignments of the coefficients c£^ , , C 3 , given by m , and define 

the polynomials p^ (z) by (flTll . The 3 zeros z^, z^\ z^ can be easily 
obtained from (1141) : their approximate numerical values are given by (1151) . Of 
course these numbers are defined up to permutations, but, as explained below, 
this has no relevance for the eigenvalues of the matrices M/ 1 ) and M*- 2 ), see 
Remark 2.2. 

For N = 3 the coefficients from (15cl) are given by csi). By inserting 
these values of the quantities into formula m for N = 3 we obtain the 










components of the corresponding 6 matrices IVf 2 )^: 


M (2 )(^) = - 


,(/d _(a») 


r„(M) m i 

Zn - Z n+1 


Z n z n +2 






Jm) 

' / m+l 


M 

^m +2 


_ J/d od 

^m +1 ~m +2 


+6 


A#d _ X 

•''M X 


/ M 

6 n 


1 4. *0d 4. z (*d 

1 r- ^ m +l T z m+2 


x i ^m +1 ^ra +2 


~(/d 


(4> _ Od 
C 1 c 2 

1 


1 4 


„(/d „(/d 


+ ■ 


J/d 

-'m +1 


» 


(aO „(/d 

^77142 "C ^m+1 ~m+2 


J#d _ (/d 

c 2 °3 



to = 1,2,3 mod(3) . (19) 


It can be checked by direct computation that the 6 matrices given by m 
all feature the 3 eigenvalues 1, 4, 9. Moreover, it can be verified, again by 
computation, that a permutation of the zeros of the polynomial (z) results in 
an appropriate permutation of the rows and the columns of the matrix IVf 2 )^, 

see Remark 2.2 . ■ 


3 Proofs 


In this Section 3 we prove Propositions 2.1 and 2.2. 

Proof of Proposition 2.1. Our starting point is the system of N ODEs 


N 


1 -1 


c„ = - 

(see eq. (9a) of [3]), with 


n (c - q ) 

= 1 , 


N 

57 [^m (CJ 


N—m 


(20a) 


7m i 


N 


7m - 57 ~ llY 

£—1, £z£m 


(20b) 


Above and hereafter we assume the 2iV quantities ( n = ( n ( t ) and 7 m = (t) 

to be functions of the independent variable t, and we denote with superimposed 
dots differentiations with respect to this variable (the dependence on which is, 
for notational simplicity, not always explicitly displayed, see for instance (1201) '). 
We moreover assume GUI the quantities r ) m = y m (t) respectively C, n = C, n ( t ) 
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to be the TV coefficients respectively the TV zeros of a time-dependent monic 
polynomial ip N {(;t) of degree TV in the variable (: 


N N 

V’iv (C; t) = C N +J2 [7 m w c N ~ m ] = n [C - m 


m= 1 


(21a) 


n=1 


implying 


7m (*) 


(—!) m 


E C ai (*K. a (*)-"C. m W 

l<si<S 2 <...<s m <AT 


(21b) 


It is this connection that justifies the validity of II 20 all : see [3111]. 

It is on the other hand well known (see for instance [6] sect. 2.3.4.1) that dy¬ 
namical system (I20bl) is integrable indeed solvable and isochronous with period 
2 -7T, all its solutions featuring the remarkable property 


7m (t + 2 7r) = 7 m (t) . (22) 


Hence, for the dynamical system implied by (I20al) with (I20bl) . reading 



N 

1 N { 

N 

\ 

Cn ^ 

II (Cn - U) 


7m - (^m - ItY 1 

(U N ~ m 


1=1, l^n 

m—1 1 

£= 1, l^m 

J 


(23a) 

with 7 m = 7 m (t) expressed in terms of the TV zeros £ n = £ n (t) via (121bl) , all 
the solutions (t) are likewise periodic with period 27r, being the N zeros of a 


polynomial ip N of degree IV in £ (see (121am . which is itself periodic with 
period 2n since its coefficients are all periodic with period 2-7 t: 


C„ (t + 2 7T) = C„ (t) • (23b) 

Remark 3.1. It might be observed that the zeros ( n (t) of a time-dependent 
polynomial of degree TV in £ which is itself periodic with period T might them¬ 
selves be periodic with a period which is a (generally small; see J 5 J) integer 
multiple of T, due to the possibility that the zeros, as it were, exchange their 
roles through their time evolution. But this is certainly not the case for a time 
evolution in which each zero evolves periodically remaining infinitely close to 
its equilibrium position in an equilibrium configuration of system (12311 ; which is 
the case we consider below. ■ 

It is moreover clear that an equilibrium configuration of this dynamical sys¬ 
tem, (1 23 all . is provided by the values £ n ( t ) = z n which are the zeros of the 
polynomial (I21al) corresponding to the equilibrium configuration 7 n (t) = c n 
of the system (I20bl) . where the N quantities c n satisfy the system of algebraic 
equations 

N 

Cm- (c m - q ) -1 = 0 , to = 1,..., TV . (24a) 

£= 1 , 
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It is on the other hand well known (see for instance Appendix C of [6]; this 
finding is actually much older, see for instance [9] and references therein) that 
the TV zeros c n of the Hermite polynomial of order TV in c, 

H n (c„) = 0 , n = l,..., TV, (24b) 

satisfy this system of algebraic equations, (124al) . 

Next, let us look at the behavior of dynamical system K3l) in the infinitesimal 
vicinity of the equilibrium configuration C, n = z n , where the TV coordinates z n 
are the TV zeros of the polynomial pn (z) , see 0- To this end we set 

Cn (t) = z n + e v n (t) , 7 m (f) = c m + £ w m (t) + O (e 2 ) , (25a) 

with £ infinitesimal; and we note that, via (I24al) with (|21b[) . these two formulas, 
(125al) . imply that 


N 

lm~ (rim-ltT 1 

i=l, 


= £ 

where 

Wj (t) = (-1) J 


N 

Wm+ ^2 ( c m ~ ce)~ {w m - we) 

1 = 1 , 


+ 0(e 2 ) , (25b) 


1 E 

j j 

it) n 

\ 


9=1 r=l, r^q 

1 


(25c) 


It is then easily seen - via these formulas and the insertion of (I25al) in (B^l) - that 
the dependent variables v n = v n (t) evolve in time according to the following 
linearized version of (l23l) : 


, (t) = -i 


N 


[J (z n - ze) 

1=1, l^n 


N N , , , 

N-j r... w 0 W ~ w t w 


J 2 ( Z n) N ° Wj(t)+ J2 i , 
j =i 1 = 1 , i^j c C 


(26) 

(27) 


where of course the quantities w n ( t ) must be replaced by their expressions (125cl) . 
Clearly this system can be more compactly rewritten as follows: 

v(l) = iM«v(l) , (28) 

with the time-independent TV x TV matrix MA 1 ) defined componentwise by (| 8 |). 

Remark 3.2. It is easy to see that a switch between the two zeros z n and 
z m of the Hermite polynomial Hn(c) in the derivation of system (f27jl . see (125al) . 
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will result in v n (t) and v m (t) exchanging their roles. This, in turn, will result 
in a switch of the n-th and the m -th rows as well as the ro-th and the m-th 
columns of the matrix M* 1 ) in system (1281) . see also © and Remark 2.1. ■ 
System (1^51) implies that the iV-vector v (t) evolves as follows: 


N 

V (t) = ^ [ a ™ eX P (i t) U (m) 
772=1 


(29a) 


where respectively A m are the eigenvectors respectively the eigenvalues 
of the N x N matrix M^ 1 ) (see (|8jl). and the N parameters a m are of course 
characterized by the initial values of the iV-vectors v so that 


N 

v ( 0 ) = ^ [a m 

772=1 


(29b) 


But we know that the time evolution of the N coordinates ( n ( t ) is periodic with 
period 2 - 7 T, see (I23bl) . hence (see the hrst of the two formulas (125aD l the iV-vector 
v(f), of components v n ( t ), must also be periodic with period 27 t. Hence—see 
(I29al) -the N eigenvalues \ m of the N x N matrix must have the integer 
values m identified in Proposition 2.1, see ©. Q. E. D. 

Proof of Proposition 2.2. Here we derive an analogous result to that 
reported in Proposition 1.1 -identifying thereby a N x N matrix M^ 2 ) (z) 
featuring as its N eigenvalues the squared integers m 2 . The derivation of this 
finding—which was actually the first one we obtained, as follow-up to 0 — is 
reported below rather tersely, since it is quite analogous to the proof of Propo¬ 
sition 1.1, see above. Also the notation is generally identical to that used 
above, although the quantities used below should not be identified with those 
used above; only in some case we have appended the upper symbol ^, to em¬ 
phasize the difference of the results reported below from those of Proposition 
2.1 and its proof. 

Our starting point is the novel dynamical system characterized by the N 
Newtonian equations of motion of goldfish type (see 0 ) 


with 


N 


L= E 


£= 1 , 

N 


2 Cn 0 

Ct2 Cj£ 


II - C<) 

£=1, i^n 


N 

[7m (Cn)^ 


(30a) 


N 

7m 7m ~b 2 ^ ( (7m T^) 

t= 1, 


-3 


(30b) 


We again assume, as above, that the quantities £ n = C, n (t) respectively 7 m = 
7 m (t) are the N zeros respectively the N coefficients of a time-dependent monic 
polynomial i/j n {(;t) of degree N in the variable see (12T1) . 
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It is on the other hand well known that dynamical system (I30bl) is integrable 
indeed solvable and isochronous with period 27r, all its solutions featuring the 
remarkable property 

7m (t + 2 tt) = 7 m (t) (31) 

(see for instance (4j [6] [7]). Hence, for the dynamical system implied by (130al) 
with (l30bD . reading 


N 


Cn= E 


1 =1, t±n 


2Cn C< 

Cn Ci 


N 


n - Cl) 

= 1 , £^n 


N 


■E 


N 

7m 2 ^ ( (7m 7f) 

i= 1, ^m 



(32) 


with 7 m = 7 m (f) expressed in terms of the TV zeros C n = C, n (t) via (I21bl) . all 
the solutions C, n (t) are likewise periodic with period 27r, being the TV zeros of a 
polynomial rp N of degree TV in ( (see (121al) l. the coefficients of which are 
all periodic with period 27r: 


Cn it + 2 7T) = Cn (*) • (33) 

Indeed an equilibrium configuration of dynamical system (15^1) is provided by 
the N zeros z n of a polynomial p n (z) the N coefficients c m of which satisfy the 
set of N algebraic equations 

N 

- Cm + 2 E ( c m - c e )~ 3 = 0 ; (34) 

£=1, £^m 

since clearly the right-hand sides of the N ODEs (l32l) all vanish at equilibrium- 
i. e. when all the “velocities” ( n vanish—provided moreover j m = c m , m = 
1,..., TV with the TV parameters c m satisfying the TV algebraic equations (l34l) . 
But it is well known (see for instance mm that the TV zeros c„ of the Hermite 
polynomial Hn (c) do satisfy the set of TV algebraic equations (1341) (in addition 
to satisfying the different set of TV algebraic equations (|24all h 

Next, let us look at the behavior of dynamical system (15^1) in the infinitesimal 
vicinity of the equilibrium configuration = z n , where the TV coordinates z n 
are the TV zeros of the polynomial p n (z) , see ©■ Then via (|25a|) we get 


N 

7m d - 2 E (7m 

£=1 , £^m 

+ 0(e 2 ) , (35) 

of course with Wj (t) defined by ()25c[) . 


= — £ 


Wm + 6 


N 

E 

£=1, £^m 


(c m - C() 4 (w m - m) 
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It is then easily seen- via these formulas— that the insertion of (I25al) in (1321) 
implies that the dependent variables v n = v n (t) evolve in time according to the 
following linearized version of (15^1) : 



N 

(<) = - 

J] ( Z n - Zt) 


e=i, e^n 

N 

f N 


\ Wj (t) + 6 Y 


J= i l ■«=!, 


l 


(t) 

{cj 



(36) 


where the quantities w n ( t ) must be replaced by their expressions (125cl) (of course 
with v m ( t ) replaced by Vm (t)). 

Clearly this system can be more compactly rewritten as follows: 

v (2) (t) = -M (2) v (2) (t) , (37) 


with the time-independent TV x TV matrix defined by (fTUb . 
System (EH implies that the N -vector v (t) evolves as follows: 


N 


v<»> (o= y : 

771=1 



(38a) 


where u^ 2 ^ m ^ respectively 



are the eigenvectors respectively the eigenval¬ 


ues of the TV x TV matrix M^ 2 ), and the 2TV parameters a™ and b^m are of course 
characterized by the initial values of the two TV-vectors v and v so that 


N 

V® (0) = Y, [ a ™ u(2)(m) 

771=1 


N 

v(2) (0) = Y [€ ] U (2)(m) 
771=1 


(38b) 


But we know that the time evolution of the TV coordinates ( n (t) is periodic with 
period 27 t, see (1551) . hence (see the first of the two formulas (l25aD l the TV-vector 
v(f), of components v n ( t ), must also be periodic with period 27 t. Hence—see 
(138al) -the TV eigenvalues Aof the TV x TV matrix MA 2 ) (z) must have the 
squared-integer values m 2 , m = 1,2,..., TV. Q.E.D. 


4 Outlook 

The findings reported in this paper suggest some further developments, which 
we intend to pursue in future publications. A quite natural one is the extension 
of these results to other, more general, classes of polynomials than Hermite poly¬ 
nomials (see, for instance, the properties of the zeros of polynomials reported in 
mm)- Another avenue of extension that we plan to explore is by iterating the 
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approach used in this paper, to find properties of the zeros of rnonic polynomi¬ 
als the coefficients of which are the zeros of monic polynomials the coefficients 
of which are the zeros of, say, Hermite polynomials; with an obviously ample 
choice of how this procedure can be further iterated and also combined with 
findings involving other named polynomials than Hermite polynomials. 
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